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Abstract 



a 

Q ■ We derive, using the algebraic Bethe Ansatz, a generalized Matrix 

, O, . Product Ansatz for the asymmetric exclusion process (ASEP) on a 

one-dimensional periodic lattice. In this Matrix Product Ansatz, the 
components of the eigenvectors of the ASEP Markov matrix can be 
00 ' expressed as traces of products of non-commuting operators. We de- 

^^ ■ rive the relations between the operators involved and show that they 

T^lj- . generate a quadratic algebra. Our construction provides explicit finite 

^^ \ dimensional representations for the generators of this algebra. 
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C : 1 Introduction 

o 

The asymmetric simple exclusion process (ASEP) that plays a fundamental 
role in the theoretical studies of non-equilibrium statistical mechanics, is a 
rS \ driven lattice gas model in which particles interact by hard core exclusion, 

c^ I The ASEP was originally introduced as a building block for models of one 

dimensional transport where geometric constraints play an important role 
(e.g., hopping conductivity, motion of RNA templates and traffic flow). 
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The exclusion process is a stochastic Markovian model whose dynami- 
cal rules are encoded in an evolution (Markov) matrix. Exact results for 
the ASEP in one dimension have been derived using two complementary ap- 
proaches, the Matrix Product Ansatz and the Bethe Ansatz (for a review see 
Derrida 1998, Schiitz 2001). The Matrix Product Ansatz (MPA) (Derrida et 
al. 1993) is based on a representation of the components of the steady state 
wave function of the Markov operator in terms of a product of matrices. This 
method has been used to calculate steady state properties of the ASEP such 
as the invariant measure (Speer 1993), current fluctuations in the stationary 
state and large deviation functionals (Derrida et al. 2003). 

The ASEP is equivalent to a non-hermitian Heisenberg spin chain of the 
XXZ type and can be mapped into a two-dimensional six-vertex model at 
equilibrium : these mappings permit the use of integrable systems techniques, 
such as the Bethe Ansatz. The Bethe Ansatz provides spectral information 
about the evolution operator (Dhar 1987, Gwa and Spohn 1992; Kim 1995, 
Golinelli and Mallick 2004) which can then be used to derive large deviation 
functions (Derrida and Lebowitz 1998). 

The exact relation between these two techniques has been a matter of 
investigation for a long time (Alcaraz et al. 1994; Stinchcombe and Schiitz 
1995). In a recent work, Alcaraz and Lazo (2004) have expressed the eigen- 
vectors of integrable quantum chains (such as the anisotropic Heisenberg 
chain) as traces of products of generators of a quadratic algebra. This Ma- 
trix Product Ansatz leads to the Bethe Ansatz equations of the system. 

In this work, we solve the inverse problem: we prove that a Matrix Prod- 
uct representation involving quadratic algebraic relations between operators 
can be deduced and constructed explicitly by applying the Algebraic Bethe 
Ansatz to the ASEP. The quadratic algebra we obtain is similar to the one 
studied by Alcaraz and Lazo. However, our algebra satisfies different bound- 
ary conditions that modify drastically the properties of its representations 
and ensure the existence of finite dimensional representations. 

The outline of this work is as follows. In section El we apply the Alge- 
braic Bethe Ansatz to the totally asymmetric exclusion process. In section 
El we derive the Matrix Product Ansatz from the Algebraic Bethe Ansatz. 
In section |3l we establish the quadratic algebraic relations satisfied by the 
operators of the MPA. Our results are generalized to the partially asymmet- 
ric exclusion process in Section El Concluding remarks are presented in the 
last section. In the appendix, we derive the Bethe Ansatz equations from 
the quadratic algebra of section HI 
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2 Algebraic Bethe Ansatz for ASEP 

We consider the exclusion process on a periodic one dimensional lattice with 
L sites (sites i and L + i are identical). A lattice site cannot be occupied by 
more than one particle. To represent the state of a site, z (1 < i < L) we 
use the spin-1/2 language : a site i can be in two states that we label as | t) 
(i is occupied) and | J,) {i is empty). A configuration C is represented either 
by a vector of the type | jf . . . f) or by a binary vector 

\C) = \n,...,T,). (1) 

where r^ = 1 if the site i is occupied and Tj = otherwise. The space of all 
possible configurations is a 2^ dimensional vector space that we shall denote 
by 5. 

The system evolves with time according to the following stochastic rule: 
a particle on a site i at time t jumps, in the interval between t and t + dt, with 
probability p dt to the neighbouring site z + 1 if this site is empty {exclusion 
rule) and with probability q dt to the site i — 1 if this site is empty. The jump 
rates p and q are normalized such that p + q = 1. In the totally asymmetric 
exclusion process (TASEP), the jumps are totally biased in one direction 
(p = 1 and g = 0). For sake of simplicity, we shall discuss the TASEP case 
in full detail. The general case will be considered briefiy in Section El 

We call ipt{C) the probability of a configuration C at time t. As the 
exclusion process is a continuous-time Markov process, the time evolution of 
iptiC) is determined by the master equation 

j^MC) = Y.M{C,C')MC'), (2) 

o 

where the element M{C, C) is the transition rate from configuration C to C 
and the diagonal term M(C, C) = — X^c M{C', C) represents the exit rate 
from configuration C. The Markov matrix M, that encodes the dynamics of 
the exclusion process, is a square matrix of size 2^ acting on the configuration 
space S of the TASEP. This Markov matrix can be expressed as a sum of 
local operators that update the bond {i,i + 1) : 

L 

M = ^M,,,+i, (3) 

where the TASEP local update operator Mj_j+i is given by 



Mi,i+i = li ® I2 . . . li-i 
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The matrix Ij is the 2x2 identity matrix acting on the site j and the 4x4 
matrix appearing in this equation is written in the local basis (| IT )? I Ti )? 
lit )i I ii ) )! that represents the four possible states of the bond {i,i + l). 
Thus, Mj^j+i is a 2^ matrix that acts trivially on all sites other than i and 
i + 1 and makes a particle jump from site i to the site i + 1. 

It is also useful to define the permutation operator Pj^j+i that exchanges 
the states of sites i and i + 1: 



Pi,i+i = li ® I2 • • • li-i 



/ 1 N^ 

10 

10 

V 1 y 



'-j+2 



(5) 



More generally, we can define a jump operator Mij and a permutation 
operator Pij, between sites i and j that act trivially on all sites other than 
i and j; Mij makes a particle jump from site i to site j and Pij exchanges 
the states of sites i and j. 

In the Algebraic Bethe Ansatz method (see e.g., Nepomechie, 1999 for 
an introduction to this subject), an auxiliary site a is introduced, which can 
be in two states labelled as |1) (site a is occupied) and |2) (site a is empty). 
These two states span a two dimensional vector space. A, the auxiliary space. 
We now define an operator £j(A) that acts on the tensor space A^ S; this 
operator acts trivially on all sites other than a and i, and is a function of a 
spectral parameter A 

A(A) = P,,,(1 + AM,,,) . (6) 

This operator can also be represented as a 2 x 2 operator on the vector space 

'^'(^) - ( :!^)' :r.! ) ' (^) 

where the matrix elements a(A), 6(A), c(A) and d{X) are themselves 2^ x 2^ 
operators that act on the configuration space S. These operators act trivially 
on all sites different from i and are given by 




(9) 
(10) 

(11) 
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The tensor products in these equations represent products over the two- 
dimensional local configuration space of a site. 

The operators Ci{X) satisfy a Yang-Baxter type relation. We consider the 
operator 

7^(z/) = 1 + uMa',a , (12) 

that acts on A^ A' where the auxiliary spaces A and A' correspond to the 
auxihary sites a and a'. In the basis (|la, !„'), jla, 2a/), |2a, !„'), |2a, 2^/)) of 
A<^ A', the operator 7l{i') is represented by a 4 x 4 scalar matrix : 



The following identity is then satisfied 



f 1 0\ 

1 z/ 

1-z/ 

V 1 / 



(13) 



7^(z/)[A(A) ® £',(/!)] = [A(/i) ® /:^(A)]7^(^.) with z/ = ^-^, (14) 

where Ci and C'i are interpreted as 2 x 2 matrices acting on A and A', 
respectively, with matrix elements that are themselves operators on S. Their 
tensor product is thus a 4 x 4 matrix, acting on A® A' with matrix elements 
that are operators on S. 

The monodromy matrix is defined as 

T(A) = £i(A)£2(A).../:l(A), (15) 

where the product of the £j's has to be understood as a product of 2 x 
2 matrices acting on A with non-commutative elements. The monodromy 
matrix T(A) can thus be written as 

where A, B, C and D are operators on the configuration space S. Taking 
the trace of the monodromy matrix over the auxiliary space A, we obtain a 
one-parameter family of transfer matrices acting on S 

t(A)=Tr^(T(A)) = A(A) + D(A). (17) 

Equation ()14|) implies that the operators t(A) form a family of commuting 
operators (see e.g., Nepomechie, 1999). In particular, this family contains the 
translation operator T = t(0) (that shifts all the particles simultaneously one 
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site forward) and the Markov matrix M = t'(0)/t(0). Using Algebraic Bethe 
Ansatz, the common eigenvectors of this family are explicitly constructed by 
actions of the B operators on the reference state Q, defined as 

n = ITT •••!)• (18) 

The state Q corresponds to a configuration where all the sites are occupied. 
More precisely, for any n < L, we define the vector 

\ZI, Z2... Zn) = B{Zn) . . . B{z2)B{zi) fi , (19) 

where zi^ Z2 ■ ■ ■ Zn are complex numbers. Because each operator B creates a 
hole in the system, the state \zi,Z2... Zn) is a linear combination of config- 
urations with exactly n holes. This vector is an eigenvector of the operator 
t(A) (for all values of A) and in particular of the Markov matrix M, provided 
the pseudo-moments zi, ^2, • • •, Zn satisfy the Bethe equations : 



Zl 



;_!)"-! J] l_ii for l = l...n. (20) 



The corresponding eigenvalue of t(A) is given by 

_ (i-A)'' + A^n?.i(^i-i) ,„, , 

Using the Bethe equations ()2U|). we find that E{X) is a polynomial in A of 
degree L — n. 



3 Derivation of the Matrix Product Repre- 
sentation from Bethe Ansatz 

In the previous section, we have constructed the eigenvectors of the Markov 
matrix by using the the Algebraic Bethe Ansatz, see equation (fTTHl . In this 
section we show that the Algebraic Bethe Ansatz also permits us to express 
the components of an eigenvector as a Matrix Product (Mallick 1996). From 
equation ()16|) . we remark that 

5(A) = (1|T(A)|2). (22) 

Using this relation, the eigenvector given in equation ()19j) can be written as 

|zi,Z2...^n) = (l|T(^„)|2)...(l|T(z2)|2)(l|T(zi)|2)f^ 

= (1, 1, ... , l\T{zn) ...® T{z2) ® T{z,)\2, 2, . . . , 2) fi 
= Tr(Q„T(z„)®...®T(2i)f]), (23) 
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where the tensor-products act on the space ^®" with n auxihary sites. The 
boundary operator Q^ is given by 

Q„ = |2,2,...,2)(l,l,...,l|. (24) 

Using the definition ()lfi|l of the T matrix, we rewrite equation (j^!^ as follows 

L 

\Zi,Z2...,Zn) = TT(QnY[Ci{zi,...,Zn)^) , (25) 

i=l 

with Ci{zi, . . .,Zn) = Ci{Zn) ® . . . ® Ci{zi) . (26) 

The operator Ci{zi, . . . , Zn) is a 2" x 2" matrix acting on A®^. The matrix 
elements of £4(2:1, ... , Zn) are operators on the configuration space S that act 
trivially on all sites except the site i. We now define the two operators D„ 
and En by the following relation 

Ci{zi, . . . , 2;„) I t) = Dn{zi, . . . , z„) I t) + En{zi, . . . , z„) I I) , (27) 

where | ]) and | [) represent the states of the site i. (For sake of simplicity, 
we are writing | ]) and | [) instead of | ]i) and | |j)). 
Equivalently, these two operators are given by 

Dn{Zi,...,Zn) = {1 \Ci{Zi, . . . , Zn) \ ^) , (28) 

En{Zi, . . . , Zn) = {l\Ci{zi,...,Zn)\^) . (29) 

The operators Dn, En and Qn are 2^ x 2" matrices acting on A'^'^ with scalar 
elements. We shall now prove some recursion relations satisfied by Dn, En 
and Qn- For ra = 1, we have 

C,h)=(j l).M.O={l '7') audQ.= (; I). (30) 
From equations (PUJI and (PHj) we obtain: 

Dn+liZi, . . . , Zn+i) = (t \Ci{Zn+l) Ci{Zi, . . . , Zn)\ T) • (31) 

Using the following identities (deduced from equations (jBHUJ) 

(T|a(z„+i) = (Tl , (T|&(^n+i) = 0, (32) 

(T |c(^„+l) = (i I , {Ud{Zn+l) = Zn+l{U , (33) 

we derive the following recursion relation : 

Dn+l{zi, . . . , Zn+l) = (q ^ j^^n+f^QJO^n (34) 



DniZi, . . . , Zn) 

En{Zl, . . . , Zn) Zn+lDn{Zi, . 



(35) 
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Thus, -Dn+i is a 2""^^ x 2"+^ matrix written as a 2 x 2 matrix built of blocks 
of size 2" X 2". Similarly, we have 

En+M,...,Zn+l) = (I ^~0^^') ^^-+{1 ?)®^" (36) 

\ En{zi, . . . , Zn) J 

These recursion relations, together with equation (jHUj) . allow to calculate 
Dn and £"„ for all n. In particular, we remark that Dn is a lower triangular 
matrix and its 2"- eigenvalues are given by 11^=1 ^l'', with e^ = or 1. 

Finally, using equation (|2l|) . we obtain 



^-^' - ( Q„ ] • ^^^) 

We now prove that the operators Dn, En and Qn provide a Matrix Prod- 
uct Representation for the components of the eigenvector \zi, . . . , Zn) of the 
Markov matrix. Using equation (^, the components of \zi, . . . ,Zn) on a 
configuration C of the system can be written as, 

L 

{C\z^, ...,Zn) = Ti[Qn liinDn + (1 - r,)E„)) , (39) 

i=l 

where Tj = 1 if z is occupied in configuration C and Tj = otherwise. Hence, 
a particle is represented by the matrix Dn and a hole is represented by the 
matrix En- 

Because of the conservation of the number of particles and holes, the right 
hand side of equation ()39|) vanishes when the configuration C does not have 
exactly n holes. If we call Xi,X2, ■ ■ ■ ,Xn the positions of the n holes in C, 
equation ^^ can be rewritten as 

{x,,...,Xn\zi,...,Zn)= TliQnDl^'^ EnD^^-^^'^ En . . . D^--^^^-^ EnD^n"""^) 

(40) 
This expression gives a Matrix Product Representation for any eigenvector 
of the Markov matrix. This expression generalizes the steady-state Matrix 
Product introduced by Derrida et al. (1993). In the Appendix, we show that 
the expression ()4U|) can be recast in the familiar coordinate Bethe Ansatz 
form. 



Golinelli, Mallick — Matrix Product Representation for ASEP 9 

4 Identification of the Quadratic Algebra 

In the previous section, we have derived the Matrix Product Representation 
from the Algebraic Bethe Ansatz by constructing exphcitly the operators Dn 
and En. We now prove that these operators satisfy some simple algebraic 
relations; more precisely, the operator En can be decomposed as a sum of n 
operators E^^^ 

n 

^n = E^^^ (41) 

that obey the quadratic relations 

E«D„ = z,DnE^\ (42) 

for ty^j : (l - 1) E«E(^) = - (l - ^) E^'^E^:^ , (43) 

and ^i^^^i*^ = 0. (44) 



In equations ()42H44|1 . the scalars Zi are arbitrary complex numbers and do 
not have to be solutions of the Bethe equations. Such a quadratic algebra 
was postulated by Alcaraz and Lazo (2004) as an Ansatz to diagonalize the 
Hamiltonian of quantum spin chains. We show here that, for the ASEP, 
this quadratic algebra can be rigorously deduced from the Algebraic Bethe 
Ansatz. Our construction is explicit and provides finite dimensional repre- 
sentations of the abstract quadratic algebra defined by relations (J12H111)- 

We shall prove the relations (PTVI44|1 by induction on n, by diagonalizing 
the matrix Dn and calculating En in the new basis. We recall that all the 
matrices with index n are of size 2" x 2^. 

For n = 1, the matrix Di is diagonal, and Di and Ei satisfy the rela- 
tions (jUmil). Now, we suppose, by the recursion hypothesis, that we have 
already diagonalized the matrix Dn, i.e., we have found an invertible matrix 
Rn such that 

R-'DnRn = K , (45) 

where A„ is a diagonal matrix with diagonal 

diag(A„) = (1, zi, Z2, Z2Z1, ^3, 2:32:1, 2:32:2, 2;32;22;i, . . . , 2„ . . . 21) . (46) 

(For n = 1, we have Ai = Di and -Ri = 1). In the new basis, the matrix En 
becomes 

^n = Rn EnRn ■ (47) 

We suppose, again by the recursion hypothesis, that we have found a decom- 
position of £n 

n 
^n = E^^^ (48) 
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such that the relations ()42fB^ are satisfied between A„ and the S^^^s. Then, 
by a change of basis, the same relations are also satisfied by D^ = Rn^nR'^ 
and E« = RnSH^R-^. 

We now show that a similar decomposition can be found at the level n + 1. 
We first construct the matrix Rn+i that transforms -Dn+i into the diagonal 
form PUI) . Using equation ((23), we take Rn+i to be of the form 



i-n ^ ^71 -1 i-n 



where An is an unknown matrix to be determined. From equations ()H5| |33l 
EH and ED), we obtain 

Rn+lDn+lRn+l = A\^-y'^\A^P ■y A ' ^^^^ 

This matrix is diagonal if and only if A^ satisfies the relation 

Anl^n ~ ^n+l^nAn = c.„ . (51) 

Knowing that A„ and £^^ satisfy the relations ()42ffl^ . we find the solution 
An of equation ^T\ : 

n c{i) 

An = K'E—^- (52) 

We thus obtain 

A„+i = Rn+iDn+lRn+l ^ \ Q ^ A ) ' ^^'^^ 

/S.n+1 is a diagonal matrix. 

The operator En+i in the new basis is found to be : 



£-n+l — Rn+l^n+lRn+l 



(1 ~ Zn+l) y , — (1 — Z„+i)A„ 



(54) 
This equation is derived by using equation (j3^ and the relations (J12H11|)- 
We emphasize that Sn+i is a strictly upper-triangular matrix : its lower-left 
elements and its diagonal vanish identically. From equation (J54j) . we deduce 
the decomposition of Sn+i 

n+l 

£-n+l = 2^ ^n+1 ) (55) 

i=l 
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where 



for i<n, Sl^l^ 



cLilCl ^ri+l 



i-^„+i)£:« 




^n+l \Zi 




Knowing that A„ and 8^"^ satisfy the relations (J12H11|), and using the 
exphcit expressions (|^ I3U1 and ISTj) . we find that the operators A„+i and 
En+i also satisfy the algebraic rules (jI^Hiij) for 1 < z < ra + 1. Reverting 
to the original basis by using the matrix Rn+i, we conclude that Dn+i and 
En+i satisfy the same relations. We have thus shown the existence of the 
quadratic algebraic relations P^HHj) at the level n + 1. 

We finally discuss the existence of relations between the operators D„ 
and Eni and the boundary operator Qn- Using the equations (j^ and (jH^j) . 
we obtain the following relation between D„ and Qn '■ 



DnQn =[l\^i] QnDn =[T\zi\Qn. (58) 





We emphasize that in the quadratic algebra that we have derived from the 
algebraic Bethe Ansatz, there is no algebraic relation between E^Qn and 
QnEn- Therefore, the quadratic algebra that we have constructed is akin 
to but not identical to that studied by Alcaraz and Lazo (2004). However, 
any modification of the boundary relations alters the properties of the alge- 
bra and profoundly modifies its representation theory. For example, it can 
be proved (Golinelli and Mallick, in preparation) that, for the ASEP, the 
algebra defined in (Alcaraz and Lazo 2004) is such that all its finite dimen- 
sional representations have vanishing traces and, therefore, can not be used 
to construct a Matrix Product Ansatz. In contrast, the algebra we have con- 
structed here admits finite dimensional representations with non-zero trace 
and therefore allows us to define a bona-fide Matrix Product Ansatz. Besides, 
from a physical point of view, it is well known that boundary conditions play 
a crucial role in the ASEP (see e.g., Schiitz 2001). 

The algebra (J12H111) and the boundary equation (jKHjl encode the Bethe 
Ansatz. In the Appendix, we use these relations to prove ab initio that 
the vector \zi, . . . ,Zn) whose components are given in equation pn|) is an 
eigenvector of the Markov matrix M provided the pseudo- moments zi, Z2, 
. . ., Zn satisfy the Bethe equations ((201) • 
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5 Generalization 

In this section, we briefly indicate how to generahze our results to the par- 
tially asymmetric exclusion process in which the particles hop to the right 
and to the left with jump rates given by p and q respectively. For the par- 
tially asymmetric exclusion process the local update operator Mj j+i is given 
by (we omit the identity operators for sake of clarity) : 



M 



i,i+l 



f N^ 

—p q 

p —q 

V J 



(59) 



Using this local operator, the construction of the £j matrices and the mon- 
odromy operator T is identical to that explained above. The recursion re- 
lations (J35|) and (J37j) for the Dn and En operators, respectively, are now 
replaced by : 



Dn+l{Zl, . . . , Zn+l) 

and 

En+l{Zl, . . . , Zn+l) 



Dn{zi,...,Zn) \ 

1 • • • 5 Zn) J 

(60) 



,1 — -Zn+l)En{Zl, . . . , Zn) Zn+lDn{Zi, . . . , Z. 



I -Zn+lEn{Zl, ■ ■ ■ , Zn) (1 — Z„+i)-D„(Zi, . . . , 2;„) \ 

V En{Zi,...,Zn) J 

(61) 

The operators Di and Ei are identical to those defined in equation (pH) . 

Here also. En can be written as a sum of n operators E^^^ as in equa- 
tion ()4H) . The operators E^^^ and Dn generate a quadratic algebra. The 
relations (jl^ and (jl^ still hold good but the equation (jIH)) is replaced by 



ioTZ^j : (l - ^ - Q^?) E^Elp = - (^1 - ^ - qz^ E^^ E^ . (62) 

These algebraic relation are obtained, again, by recursion on n. The diag- 
onal basis is found by using the transformation i?„ defined recursively in 
equation P^ where the matrix An is now given by 

An = {l- ^Zn+,)K' E -^ • (63) 

P j^^i Zi — Zn+l 

For g = 0, we recover the expressions given in sections El and IH 
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6 Conclusion 

In this work, we have shown that the components of the eigenvectors of the 
asymmetric exclusion process can be written as traces over matrix products. 
This Matrix Product Representation has been constructed from the Alge- 
braic Bethe Ansatz in a systematic manner ()40j) . Our method also allows to 
derive the algebraic relations ()41fH^ satisfied by the operators that repre- 
sent particles and holes. The quadratic relations obtained are in fact logical 
consequences of the Algebraic Bethe Ansatz procedure and thus, ultimately, 
stem from the Yang-Baxter equation. The approach described in this work 
shows that there is a close relation between the Matrix method and the Bethe 
Ansatz, at least in the case of the ASEP. We believe that the derivation of 
the Matrix Ansatz for the ASEP presented here can be generalized to inte- 
grable quantum chains. Besides, if the equivalence between Matrix Ansatz 
and Bethe Ansatz is true in general, this would provide a technique for con- 
structing the quadratic algebras from first principles rather that having to 
postulate them a priori. We also emphasize that, in contrast with the work 
of Alcaraz and Lazo (2004), our construction provides explicit finite dimen- 
sional representations of the algebras involved that can be used for actual 
calculations on finite size systems. 
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Appendix: Derivation of the Bethe Equations 

We explain here how to derive the Bethe equations (j^UI) from the quadratic 
algebra ()42fB^ and the boundary relation (j3Hj) . We first show that the 
expression (pn|) is equivalent to the standard coordinate Bethe Ansatz form: 

\Xi, . . . , Xn\Zi, . . . , ZriiJ 

where a belongs to S„ the permutation group of n objects. This formula is 
obtained by inserting the decomposition (|^ for the operators E^ and by 
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noticing from relation ()44|) that each E^^ must appear only once. We use 
equation (jl^ to push all the operators E^^'^''^^ to the right. We thus obtain 

\X\^ ■ ■ ■ 1 Xfi\Zi, • • • 5 Zjif 

E ^'[uT^ln-iT" ■ ■ ■ 4ir''"Tr(g„D^"E('^(i))E('^(2)) . . . Ei'^W)) . (65) 

We use equation (jlH)) to rearrange the product E^^^^ . . . -E^^"-* in the canonical 
order E!^^ . . . E!^^ and obtain 

{x„ ...,xr,\z„...,z^)=KY: i-ir n {^i) - ^y ^.Ti) • • • ^;w (66) 

(7^'Sn 2 = 1 

with 

K= Il^^] n(l--)"'Tr(g„Dr"E«Ef ...eW), (67) 



\i=i J i=\ 



Zi 



and where (—1)°^ represents the signature of the permutation a. We thus 
find that the eigenvector can be written (Golinelli and Mallick 2005) as a 
determinant of a matrix Vij : 

(xi, . . . , Xn\zi, ...,Zn)=K det(yij) with Vij = {zj - lyzj''' . (68) 

We now show that for the vector ()4()|1 to be an eigenvector of the Markov 
matrix M, the pseudo-moments zi, . . . Zn must satisfy the Bethe equations pUjl . 
Using the fact that an eigenvector of M is also an eigenvector of the trans- 
fer matrix t(A) for any value of A and therefore of the translation operator 
T = t(0), we obtain 

[Xi, . . . ,Xn\T\Zi,Z2 ■ ■ ■ Zn) = {Xi — 1, . . . , Xn — ^Zi, Z2 . . . Zn) 

= ({xi,...,Xn\zi,Z2...Zn) , (69) 

with ^-^ = 1 . We now substitute equation pUjl in this identity. We have to 
distinguish two cases: a;i > 1 and xi = 1. For xi > 1, we have 

Tl(QnDl^''^EnDl^"'^~^En . . . D^"-^"-^~^E„£'^^"^'"D 
C TT{QnD:,'~'EnD:'~^^-'En . . . D^,"-^--''' EnD^~^" 

(f[ z] 1i(QnDl--'EnDl^----^En . . . D:"-^"--iE„D^^") . (70) 

To derive the last equality, we have used equation (jHHl)- We thus have 

Zi...Zn = C,- (71) 
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For xi = 1, equation (j69|) becomes 

Tr(Q„D^^-2E„ . . . D^"--"-i-iE„D^--"E„ 
C Tr (QnEnD^^'-^E^ . . . ^^^--^-^-^E^D^^") . (72) 

Using the decomposition ()(i4|l . we obtain the following sufficient condition for 
equation (f?^ to be satisfied by any a : 

Tr(n r)^'2-2 E^(o-{2)) p)Xn-x„-i-l rp pfL-Xn pMl))\ _ 

J.L \^Vn-^n -^n • • • ^n ^n-^n ^n ) ~ 

Using equation P3|l . we commute the operator E)^^^'' with all the other 
operators E^^^^\ for j ^ 1, and bring it back to the rightmost position. This 
leads to the consistency condition 

From relation (f7T|l . we conclude that this equation is identical to the Bethe 
equation (j^U)). 
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